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Abstract 

Solutions to scalar theories with derivative self-couplings often have regions 
where non-linearities are important. Given a classical source, there is usually a 
region, demarcated by the Vainshtein radius, inside of which the classical non- 
linearities are dominant, while quantum effects are still negligible. If perturbation 
theory is used to find such solutions, the expansion generally breaks down as the 
Vainshtein radius is approached from the outside. Here we show that it is possible, 
by integrating in certain auxiliary fields, to reformulate these theories in such a way 
that non-linearities become small inside the Vainshtein radius, and large outside it. 
This provides a complementary, or classically dual, description of the same theory 
- one in which non-perturbative regions become accessible perturbatively. We con- 
sider a few examples of classical solutions with various symmetries, and find that in 
all the cases the dual formulation makes it rather simple to study regimes in which 
the original perturbation theory fails to work. As an illustration, we reproduce by 
perturbative calculations some of the already known non-perturbative results, for 
a point-like source, cosmic string, and domain wall, and derive a new one. The 
dual formulation may be useful for developing the PPN formalism in the theories 
of modified gravity that give rise to such scalar theories. 



1 Introduction and summary 



Perturbation theory is often the only analytic tool available to extract detailed information 
from interacting theories. The regime in which perturbation theory is valid is usually 
limited. In certain cases, however, it is possible to reformulate a theory in terms of 
new, dual, variables that allow perturbative calculations in the regime where the original 
formulation was non-perturbative. 

In this note, we discuss certain special nonlinear theories, and show that at the classical 
level they admit a dual description. These are field theories of a scalar, </>, with purely 
derivative nonlinear terms, that nevertheless give equations of motion with no more than 
two time derivatives. 

Our main motivation for considering such models, and their classical duals, stems from 
the theories that modify General Relativity (GR) in the infrared - the five- dimensional 
DGP model [1], and four- dimensional ghost-free massive gravity [21 13]. The four-dimensional 
scalar Lagrangians discussed here capture parts of the full gravitational theory, as shown 
in [1] for DGP and [2] for massive gravity (for reviews and experimental limits, see Refs. 
[7], and for a recent theoretical review of massive gravity see Ref. [S]). 

Our analysis may have broader applications though: derivatively self-coupled theories, 
in particular the galileons [H El |6], can also be obtained in the probe-brane limit of higher 
dimensional constructions [9], and their extensions [101 El [121 [ISl Ej, while their three- 
dimensional counterparts are obtained in the context of three-dimensional "new massive 
gravity" [13], as shown in |16] . 

In the present work we shall focus on the so-called cubic galileon, ~ □0((90)^/A^, [HE] 
by which the free Lagrangian, —(90)^/2, is supplemented in four-dimensions. The state 
described by can be thought of as a Nambu-Goldstone boson, nonlinearly realizing (a 
limit of) broken higher dimensional Poincare or diffeomorphism invar iance 

C + C^X^, 

with denoting a constant vector. In parallel, we will also consider - mainly as toy exam- 
ples - theories of an "ordinary" Nambu-Goldstone (NG) boson with the self-interaction 
terms, such as — (50)^"/A^"~'^, n>2. 

The only physical scale in these models is A. Such theories are usually regarded as 
effective field theories valid at energies/ moment a below A. One reason for this is that 
scatterings of the 0-quanta - when treated in conventional perturbation theory - exhibit 
non-perturbative behavior at/above the scale A. On the other hand, the galileons do not 
seem to represent garden variety effective field theories. They are special - for instance, 
they do not get renormalized by quantum corrections [H O [6l [10] (although, other higher- 
derivative terms may be generated). One may wonder then, if there may be some hidden 
structure in the galileon theories that would enable one to deal in a controllable way with 
scales above A, by a re-summation of perturbative diagrams, or by a dual description. 
Although, in the present work, we will not explore the above important question, we'll 
make a step in that direction. 
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What we shall show, instead, is that a dual description is possible for these theories in 
the classical regime. To enable a window in which the classical description is meaningful, 
we introduce a minimal coupling of the scalar cj) to the trace of an external classical 
stress-tensor T (planets, stars, etc.) 

^ T, (1) 



Mpi 

where Mpi is the Planck mass. Coupling to the stress-tensor can be non-minimal in a 
more general context [2], with certain interesting observational consequences [17]; our 
analysis should straightforwardly apply to those cases too. 

The presence of the Planck mass, in addition to A, gives rise to a new derived scale, 
referred to as the Vainshtein scale [TS]. For a static spherically symmetric classical source 
of mass M, ~ A~^(M/Mpi)^/^ ^M- This scale is much greater than A~^. The conven- 
tional perturbative expansion can be used to compute the field configuration outside the 
Vainshtein radius, r >> r^,. Inside the radius, r << r*, classical non-linearities in are 
dominant, and the perturbative expansion breaks down. More formally speaking, external 
classical sources introduce a new expansion parameter, that captures the strength of 
classical nonlinearities; for the galileon, ad = d'^(f)/A^, while for the NG-type theories, 
ad = (50)^/A^. The parameter ad is source-dependent, and for theories considered here, 
there is generically a broad region in space where > 1, while energies and momenta are 
still well-below A. Therefore, the classical field enters a highly nonlinear regime, while the 
quantum corrections are still negligible, as long as we stay at the distance scales greater 
than A"^ 

We show how these theories can be dualized by integrating in certain auxiliary vari- 
ables. The dual theory is classically equivalent to the original one, however, it no longer 
has any higher-dimensional derivatively-coupled terms. Instead, the dual theory is non- 
local, in a sense that it contains lower dimensional non- derivative terms with fractional 
classical dimensions. 

Perturbation theory in the dual version has a regime of validity opposite to the original 
one: there is still a Vainshtein radius, but now non-linearities are small inside the Vain- 
shtein radius, and large outside of it. Hence, the non-perturbative regime in the original 
variables is perturbative in the dual picture, and vice versa. 

We point out that, whether in the dual description or the original description, in both 
the perturbative and non-perturbative regimes the classical fields are weak in Planckian 
units, (pd << Mp\. This should be compared to the case of a black hole in GR, where 
nonlinearities near the Schwarzschild radius, r^, are due to classical fields that aren't 
small in Planckian units. These nonlinearities can be re-summed into the Schwarzschild 
solution^. 

We would like to make a few important comments on the literature. First, our dual 
description is what captures the properties of the small-mass expansion used by Vainshtein 
[18] in massive gravity. The latter expansion is replaced here by a series governed by 



^It may be interesting to attempt to dualize GR along the lines discussed here. 



2 



positive powers of A. This parameter, in the context of massive gravity, is derived from 
the graviton mass and Planck scale. Second, the derivatively coupled theory of a 
field containing a nonlinear ghost (for instance the Lagrangian ~ (□(/))^/A| describing 
the longitudinal mode of earlier, ghostly massive gravity theories), was decomposed by 
Deffayet and Rombouts [20] by means of the Ostrogradskii method to manifestly exhibit 
the ghost in the linear theory. Our construction is similar, but not identical, since we're 
dealing with the theories without ghosts; this essential distinction gives rise to significant 
differences in the two will be seen below. In spite of the differences, we have 

been inspired by both [18] and [20] . 

The work is organized as follows: In Section 2 we consider the simplest NG-type 
model. We discuss its classical dual and perform calculations for sources with spherical 
and cylindrical symmetry. This serves demonstrational purposes, as the main focus of 
the present work is on galileons which are relevant to the theories of IR modified gravity. 
In Section 3 we study the cubic galileon. Again, we present the dual theory, and use it 
for perturbative calculations with spherical, cylindrical, and planar symmetries. Section 
4 contains conclusions and outlook. In the Appendix we discuss more general NG-type 
theories. 

2 Nambu-Goldstone type theories 

As mentioned before, the galileons [11 16], are perhaps the most remarkable derivatively 
self-coupled scalar field theories: their special structure guarantees a good Cauchy for- 
mulation, as well as non-renormalization of these terms [H Ej in the quantum theory 
with no additional fields. Quantum effects may generate other derivative terms, such as 
(□0)*^, k > 2, however, the effects of the latter are suppressed in the classical regime 
considered here. Therefore, our restriction to a single cubic galileon term in the next 
section can be justified even in the full quantum theory. 

On the other hand, there is no similar argument for the NG-like theories. Also, there 
is no known principle which would lead to the re-summation of the ((90)^" - type NG 
interactions, except in the case leading to the DBI action. DBI, however, has no well- 
behaved static solutions [25], and we will not consider it here. Therefore, restricting only 
to the NG-type term with = 2, as it will be done in this section, is not justified in 
the full quantum theory. Nevertheless, we consider this example as a starting point in 
this section and regard it as a toy model where calculations are easier, keeping in mind 
that generically one should be retaining terms with all possible integer values for n (other 
values of n are considered in the Appendix). 

Thus, we consider a theory of a scalar field with a NG-type derivative quartic self- 
interaction 

Around the trivial background, T = 0, = 0, perturbation theory for the amplitudes 
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of the 0-quanta starts to fail as energies reach the scale A. We now turn to nontrivial 
classical backgrounds with nonzero T q 
The equation of motion reads 

o<i^+Y^d,mrdy] = -^- (4) 

For a static point-like source, T = — M(5^(x), the solution for (f) is spherically symmetric 
and static, and the equation of motion reduces to 

V0+^(V0)^V0) =^5'(x), (5) 

which can readily be integrated once to obtain a cubic algebraic equation for the radial 
derivative 0', 

This can be solved exactly. The exact solution has two regimes, depending on which of 
the two terms on the left hand side of ([6]) dominates. The scale that separates the two 
regimes is denoted r^,, this being the distance at which the two terms become comparable. 

At scales larger than r^,, the linear term on the l.h.s. dominates, leading to the usual 
Newtonian potential for the scalar 

M 1 

— , r ^ r^,. (8) 

^ 47rMpr' ^ ^ 

Note that the value of the classical field is small in Planckian units: ttf- ~ — << 1, 
for r > r*. At distances shorter than r* on the other hand, the non-linear term is more 



^Note that we consider this model with a "wrong sign" in front of the nonUnear term - the sign that 
does not admit a conventional UV completion [21], while it exhibits the Vainshtein mechanism [22]. We 
thank Lasha Bcrezhiani for pointing out that the "right-sign" nonlinear term for the NG, +{d(j>)'*'/A'^, 
does not admit the Vainshtein mechanism. 

The theory ([3]) with T = can be regarded as the decoupling limit of a massive abelian vector field 
with a quartic interaction [23], 

C = -\f^'"F^^ lm\A^A, ^ {A^A.f . (3) 

At energies parametrically above niA, 4> describes the helicity-0 component of extracted through the 
Stiikelberg replacement, A^ — > A^ — ■:^d^(j), and the decoupling limit is defined as follows, tua — > 
0, g — > 0, A = ^ fixed. The effective theory ([2]) is valid at distance scales A^^ ^ r <C mj^^. 
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important, and the solution reads 



j A^'/V^/^ + const., r<n. (9) 



It is straightforward to check that even in this regime the value of the classical field 
is sub-planckian, ~ ^{-^Y^'^ << 1; ^^(^ decreases with decreasing r. 

If the exact solution were not known, and we wanted to set up a perturbation theory 
to find it, we could perform an expansion in powers of the interaction. We could do this 
by expanding the field into powers of the non-linear interaction, 

= 00 + 01 + 02 + ■ ■ ■ . (10) 

Plugging this expansion into the equation of motion (jlj) and equating powers of j^, one 
generates a series of equations 

□00 - (11) 

+ X^^^ [(^0o)'^'^0o] = 0, (12) 
: (13) 

For the static point-like source T = — M(5^(x), the leading order solution is 0o = 

By simple power counting one can see that the series is nothing but an expansion in 

powers of the parameter 



(14) 

r / 

where is the Vainshtein radius ([7]). The expansion is good for large radii, and starts to 
fail as we approach the Vainshtein radius from the outside. 

2.1 Dual formulation 

The theory in the form given above naturally yields a perturbative expansion which is valid 
in the IR but breaks down in the UV. As discussed above, the region of transition between 
the classically perturbative and non-perturbative regimes lies around the Vainshtein radius 
I. 

We would like to rewrite the theory in a form that makes it possible to perform a 
perturbative expansion which is valid in the UV, rather than in the IR. It is straightforward 
to check that the following Lagrangian written in terms of the original scalar 0, and an 



•^We mention again that quantum-mechanically, we should really be considering all operators of the 
form (i9<^)^", which become of the same order as (90)^ once the Vainshtein radius is approached from 
the outside. Since we are concerned here with a mere illustration of how dual theories work (before a 
much more stable analysis of galileons in the next section), and for simplicity, we choose to ignore these 
operators. 
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auxiliary vector field t/'^, 

L = ~\{dcl>f + l^^'M^ir' - rd,c^ + ^0T, (15) 

recovers back the original theory ([2]) upon integrating out ip^. 

By simply introducing a new variable, ip^, we seem to have arrived at an equivalent 
action in which none of the coupling constants (never mind Planck's mass) are of negative 
mass dimension, naively pointing towards the possibility of a perturbative expansion which 
is good in the UV. Note however that there is now a potential term for if)^ with a fractional 
power, and that this fraction is less than two (in Appendix |Al we show that more general 
interactions also lead to such terms and that the fractional power is always less than two). 

In the limit A — )■ 0, Mpi — oo, the theory (fT5|) reduces to the one governed by the 
non-dynamical equations, d^cp = 0, 9^-?/'^ = 0. 

Note that the second term in (llSp . if regarded as a potential, has a "tachyonic" sign. 
This is a consequence of the minus sign in front of the nonlinear term in ([2]), and is a 
cause of the existence of superluminal modes [21], [6l [2ll [25] in the theory ([2|). Thus, in the 
dual version the superluminality is related to the tachyonic instability of the non-analytic 
potential. The latter could be stabilized, e.g., by supplementing f[T^ with carefully chosen 
higher powers of ip"^. However, we will not pursue this completion here, since we just use 
the theory ([Sj) as a toy example to demonstrate the trick of dualization in an easier setup. 

The equations of motion that follow from (ITSl) are 

□0 + d^ij'' = — —T, 

At this point we choose to make the field decomposition 

V^o = X, i/Ji = '^I + 911^, a = (j) + ilj, (17) 
under which the equations become 

A'/' i-x' + ii^I + d.i'fy'^' ii'f + d.ij) + d,i, = d,a. (18) 

One can see that the vanishing transverse component of ip, ipj = 0, is a consistent ansatz 
for the solution. Moreover, for static field configurations, x = ^ind a obeys a linear 
equation sourced by T, while a in turn determines ip. Instead of the irrelevant operator 
of the original theory ([2|), we have only the self-interaction of the field ip which looks like 
a relevant operator, and we could expect it to be subdominant in the UV. In the IR, on 
the other hand, things are ill-defined because the interactions above the trivial ground 
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(16) 



state are non- analytic]^. 

For the static point source T = —M6^{x.), the only excited degrees of freedom are a 
and the longitudinal component of ipi. The first equation in (1181) tells us that the exact 
value of a is its linear Newtonian value, a = while the equation for if) reduces 

to the following one 

If we ignore the non-linear interaction by setting A = 0, ^/^ will have the zeroth order 
linear solution, ifjQ = — ^^[j^ ^ ■ The full solution for the original field (po is then trivial, 
since there is a cancellation (po = a — tpo = 0. It is easy now to estimate the distance for 
which the non-linear term in ( !T9|) becomes important. We find that this scale is again the 
Vainshtein radius r^, ([7]). However, in contrast to the original formulation of the theory, 
the self-interactions are now important only at distances larger than the Vainshtein radius, 
r > r*. 

We can see this reversal in the region of strong coupling more explicitly by solving 
the ip equation perturbatively using the dual formulation. We set up the expansion by 
expanding in powers of the interaction coupling A^/^, 

= V^o + V^i + ^^2 + ■ ■ ■ , (20) 
plugging into (fT9l) and equating powers of A. The solution to lowest non-trivial order is 



M 1 .4/, / M -I/O 

*»+*' = -«jp;-3^' (s^Mp) (21) 

while a has the linear Newtonian 1/r profile to all orders. Recalling the definition of the 
physical field 0, we have, 

(P = a-^ = 3A^/3 + const. + ■■■ , (22) 

\47rMp / 

this shows that the expansion is in powers 

4/3 



(23) 



This expansion is inverse to the expansions of the original theory ( |T4l) . As a result, the 
dual expansion breaks down as we approach the Vainshtein radius from the inside. 



"'AH the above statements refer to the classical theory. We point out. however, that non-covariant 
decomposition (jl7p . which does not introduce extra time derivatives, is likely to be a good starting point 
for quantization of this theory. 
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2.2 Profile of an infinite string 



The use of the dual formulation is not restricted to spherically symmetric static solutions. 
It should provide a complementary description for any classical solution. For illustrative 
purposes, we now consider a cylindrically symmetric solution sourced by a uniform string, 
with the mass-per-unit-length denoted by k 

T = -M. (24) 
The exact solution again has two regimes, separated by a Vainshtein radius 

The leading behavior of (j) in the two regimes is, 

1/3 




A'^/^r^/^ + const., r <C r*, 

(26) 

r > r*, 

where Vg is a UV regulator scale - the transverse width of the string in this case. Using 
perturbation theory in the original formulation, we recover the logarithmic profile for 
r ^ r* as the leading term, and perturbation theory breaks down as we approach the 
Vainshtein radius from the outside. 

On the other hand, the dual form of the equations of motion, ( IT8l) . yields the following 
expressions for the fields a and ip, 

X = 0, a= ^ In (—^ , i) ~ ^ In ( - - ( ^ ^ A^/^r^/s ^ ^^nst. 
^ ' 27rMp \rsj 27rMp \rsj 2 \2tiMpJ 

(27) 

The expression for a here is exact as above, while the series for breaks down in the IR, as 
the Vainshtein scale r^, is approached from the inside. Recalling the definition (p = a — ip, 
one again finds an agreement with the result obtained in the original formulation of fl26|) . 



3 The cubic galileon 

We have illustrated how the simplest model of a single Nambu-Goldstone scalar can 
be re-written in a form for which classical perturbation theory has a region of validity 
opposite to that of the original formulation. The method however is quite general. The 
essence of the method is to introduce auxiliary fields in such a way as to replace the non- 
renormalizable derivative interactions with (generally non-analytic) non- derivative terms. 
We now consider the cubic galileon [H |5] - an example of a scalar field theory that is 
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relevant for modifications of gravity, 



0r. 



(2^ 



As mentioned before, in the Mpi — )■ oo limit, the galileon term in the above action does 
not get renormalized in the full quantum theory; also, no higher galileons [6] will be 
induced if they're not introduced to begin with. Moreover, no NG-type terms, (90)", 
of the previous section will be generated since the latter do not respect the "Galilean 
symmetry" — )■ + c + c^x'^. 

In spite of the presence of higher derivatives, the equations of motion that follow from 
the above action are second-order, leading to a well-defined Cauchy problem, and the 
absence of additional ghostly degrees of freedom. 



T 



(29) 



Concentrating again on radial profiles for sourced by a point-like source, T 
— M5^(x), the equation of motion fl29|) reduces to the following. 



V- 



V0 + ^ (2V0 - V(V0)^ 



M 



5'(x). 



(30) 



Integrating once, we obtain a quadratic equation for the radial derivative of the galileon 
field, 

0' 4 M 1 

Like in the model of the previous section, the exact solution has two regimes, separated 
by the Vainshtein radius r*. 



0(r) 




where 



M 



1/3 



A' 



+ const. 



r > r*. 



(32) 



(33) 



The Vainshtein mechanism is therefore at work in the cubic galileon theory as well, screen- 
ing the scalar potential significantly within r*. As in the previous section, the classical 
field is weak (sub-Planckian) in both the linear and non-linear regimes, (0/Mpi) << 1. 

Perturbation theory in this formulation allows us to compute the corrections to the 
1/r solution for r 3> r,, via a 1/A'^ expansion of the equation (!29|) . We thus write. 



(34) 
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so that after plugging into the equation of motion and equating powers of 1/A, we find 
the series of equations 

(35) 

(□0o)'] = 0. (36) 

(37) 

This gives an expansion in powers of (r^/r)^, which is vahd outside the Vainshtein radius 
and starts to fail as the Vainshtein radius is approached from the outside. 

3.1 The dual galileon theory 

We would now like to find a dual formulation of the galileon, one whose classical perturba- 
tive expansion is valid inside the Vainshtein radius. Starting with the original Lagrangian 
(128|1 . we introduce two auxiliary scalar fields 6^ and A and write an equivalent version of 
the theory as follows, 

C = -^{d<j>f + A^/'y^ - b^d,<P - Xn<f> + ^<f>T. (38) 

Again, we have succeeded in representing the cubic galileon in a form in which all terms 
look relevant at the expense of introducing fractional dimensions on fields. This has a more 
complicated structure than the dual Lagrangian of the previous section ( fTSj) . However, 
there are some similarities, such as that the nonanalytic potential term in ( l38l) also has 
a tachyonic sign, and the fact the field becomes trivial in both ( fT5|) and ( 138|) in the 
limit A — 0, Mpi — )■ oo. This, and the spectrum of the theory, is best seen bv using a 
non-covariant (i.e., a 3+1) decomposition of the vector field 6^, as given belo\MJ- 

The equations for (p, A and b'^ that follow from the latter Lagrangian are given as 
follows. 




Resorting again to the 3 + 1 decomposition of the vector field, taking the divergence 

^As noted in the previous section, this decomposition, unhke the covariant one, does not introduce 
artificially extra time derivatives which would have confused the counting of propagating modes already 
in the classical theory. 



□00 



T 



□01 - -Tj [{^^l^u<i)o) 
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of the last of the equations of motion ( !39|j I. forming suitable linear combinations and 
(re) defining fields as follows, 

bo = (3, bi = bj + d,b, X = X + b, = + A, (40) 

one can reduce (15^ to the following system, 

- do/3 + dob = -^T, 



M, 



p 



do I doX - A^/^/^/3 ]-dAdrx- + d.b) 



+dlb-doP = -^T, (41) 



As in the example of the previous section, in the presence of a point source, T = —M6^{x.), 
the transverse component bf as well as /3 vanish; moreover, one combination of the fields, 
denoted by 0, is free and has a linear equation everywhere in space, receiving an exact 
Newtonian profile 

^ ' (42) 



47rMp r 

The last equation from the system (HT]) on the other hand gives 

b = -X- irA', (43) 
which, when plugged into the second of (14T!) . gives an equation for A after integrating 



- A' + A'/\ --rX' = (44) 
V 4 47rMpr2 ^ ^ 

We have now arrived at an equation which achieves the goal of the dual formulation. 



^For the static profiles at hand, retaining only the divergence of this equation is sufficient for obtaining 
the complete solution. 

^One has to be careful with the square root at this point. Jumping a bit ahead, we note that (in 
complete analogy to the Goldstone-Stiickelberg case considered above) the leading term in the expression 
for A well within the Vainshtein radius is of the form A/r, with A denoting some positive constant. This 

uniquely fixes all signs in front of square roots, as well as makes the \/A + b expression in these equations 
well-defined. 
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The interaction term, proportional to A'^/^, becomes important only at distances larger 
than the Vainshtein radius We again set up perturbation theory by expanding the 

A-profile in powers of A^/^, 

A = Ao + Ai + A2 + . . . . (45) 
Plugging the expansion into (Hil) and equating powers of A, we obtain 

Ao = — — -, Ai = A-^/M— — — r^'^ + const., ■■■ 46 
Finally, recalling the definition of the physical field (f), we have 

/ \l/2 

= + A = —^A^rl ( — ) + const. + ■ ■ • , (47) 

in complete agreement with the result (!32l) of the original theory well within the Vainshtein 
radius. The perturbative expansion in the dual formulation is an expansion in the ratio 

(;)"■ 

and so the expansion is valid inside the Vainshtein radius, complementary to the expansion 
in the original formulation which is valid outside r*. 



3.2 Domain wall and infinite string 

Similarly to the case of the Nambu-Goldstone theory, the dual formulation should be useful 
in reorganizing the perturbation expansion for any classical solution. It is interesting to 
see this on examples other than that of a point source - such as a domain wall or a string. 

It has been shown that domain walls do not possess a Vainshtein scale in DGP and 
massive gravity (this scale is of the order of the wall width), so they give rise to a 
fifth force at all distances. This fact should be captured by the cubic galileon theory (1281) 
and therefore by its dualized version, presented above. The absence of an r^, distance 
can be easily seen in the equations of motion of the original theory f l29|) . Indeed, in the 
presence of an infinite domain wall in the x — y plane at z = 0, the problem becomes 
one - dimensional, with depending on a single coordinate z. One can then easily see 
that all nonlinearities in the original equations of motion vanish and one is left with 
the Newtonian profile of a ID source for the scalar. On the other hand, the last of the 
equations of motion (14T]) that follow from the dual theory implies A = for only z - 
dependent profiles. Recalling the definition of the original galileon in terms of the free 
scalar of the dual theory and A, = + A, one obtains agreement between the two 
representations of the theory (as must be the case since the representations are classically 
equivalent). Both the original and dual formulations are free of non-linearities, and so 
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there is no perturbative expansion to be done in either case. 

Next consider an infinite string source, T = —k^^. In the original gahleon theory, 
the equation of motion can be integrated to yield the following algebraic equation for the 
radial derivative of the axially symmetric (p profile. 



JS? r 27rMpr' 
One can immediately read off the Vainshtein radius from this equation [^, 

1/2 



(49) 



K 



MpA3 



(50) 



Using perturbation theory, one can readily solve for (p well outside the Vainshtein radius 
to obtain the leading behavior of the scalar profile, 

K f T \ 

In — , r>r=,, (51) 



where again is a cutoff scale the finite thickness of the string. Well within the Vainshtein 
radius, the nonlinear term in ( H9ll dominates, leading to 

1/2 

r + const., r <C r*. (52) 



V47rMpy 



In the dual theory on the other hand, we can proceed in complete analogy to the above 
analysis. The last of the system (HT!) yields the following identity for such profiles, 

h = -\- (53) 

Plugging this into the second of these equations and integrating, one finds that the equa- 
tion determining A is given as follows, 

-A' + (ApyiT^^^. (54) 

One can now solve this equation perturbatively, but perturbation theory now works well 
within the Vainshtein radius, 

A = In f'— ^ + ( — r + const., (55) 

27rMp \rj \AtiMv 



^Note that this decouphng-hmit expression for the Vainshtein radius coincides with the one derived 
in the full DGP model in [27l. 
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arriving at an expression for which is vahd in the UV, 



) 



1/2 



= + A = A3/2 



AirMp 



r + const. 



(56) 



This is in complete agreement with the expression (132]) . obtained from the original for- 
mulation. 

4 Conclusions 

We have studied a few examples of classical dualization for Nambu-Goldstone and galileon 
theories, which allows for a perturbative formulation of the regimes in which the original 
theory becomes classically non-perturbative. Quantum mechanically, such a formulation 
is only valid in special cases, where certain symmetries make it possible for classical 
nonlinearities to be strong, while keeping quantum corrections under complete control. 

Among scalar field theories, galileons perhaps represent the most remarkable examples 
of this, due to a powerful non-renormalization theorem [U [6l [10] protecting the leading 
part of the Lagrangian from quantum corrections. Hence, our results are justified for 
these theories. Moreover, it should be possible to generalize our approach to the higher 
galileon terms. 

In addition to capturing many features of the DGP model, galileons have emerged as 
an essential ingredient in the recently formulated ghost-free massive gravity models [21 [3]. 
The decoupling limit of these theories represents a certain (scalar-tensor) extension of the 
galileon with a more general structure, which however retains all the nice properties of the 
galileons, such as the presence of the Vainshtein mechanism and the non-renormalization 
theorems (for studies of cosmological and spherically symmetric solutions in the decou- 
pling limit and beyond in ghost-free massive gravity, see [281 1291 |30l [311 [321 [331 [311 [351 
[36l [371 [381 [391 [10] references therein). Our method should have a straightforward 
algorithmic generalization to those more general Lagrangians obtained in [2], and should 
be useful for the development of the analog of the Parametrized Post Newtonian formal- 
ism in massive gravity [3] (or for generic modified gravity models with extra scalar fields, 
such as the recently proposed Fab Four theory, [H]) and for systematically determining 
the observational consequences of the Vainshtein mechanism [121 [13]. 

The method presented above might as well be useful in studies of the proposal of 
Ref. [22] (see [H] for treatment of the scattering problem in this context for the theories 
considered above). 

Finally, it remains to be seen if quantization of the classical duals considered in the 
present work can lead to duals of the full quantum theory. Given the non-analytic nature 
of the dual theories that we obtained, quantization seems to be a nontrivial task. We 
expect the non-covariant decomposition of the auxiliary fields used in Sections 2 and 3 
to be a good starting point for bookkeeping of the degrees of freedom. The quantization 
procedure, may or may not force us to introduce additional dynamical degrees of freedom 
at the scale A. 
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A More General Interactions 

In this appendix we study a more general derivative self-interaction which is an arbitrary 
power of the field 0, in order to show that the dual formulation always has non-analytic 
self-interactions with fractional powers of the fields, and that this fractional power is 
always less than 2, i.e. less than that of a mass term. 

Consider a (ghost-free) Goldstone - Stiickelberg - type theory with a shift symmetry 
and a Z2 symmetry. A general interaction term containing only one derivative per field 
is given by the following, 

L = -lw)=-(M!): + ^^r. (57) 

The theory can be equivalently rewritten by integrating in the vector field 

L = -1(90)2-^^9,0+ (^1 - At^ + _1_^T. (58) 

Integrating out ipf^ recovers ( !57|) . Note that the power of ipfj_ in the interaction term is 
always < 2. 

The equations of motion that follow from ( l58l) can be reduced to a form similar to 
the system (fT9|) . Using the same decomposition of the vector field as before, -0, — >■ 
(x j'^f + one finds that x = "^f = 0. One combination of the fields, a = (p + ip, 
in the presence of a point-like external source has the usual 1/r profile at all distances, 
whereas i{j has a crossover Vainshtein scale due to nonlinearities. Moreover, the form of the 
dual action suggests that at small distances, nonlinearities in ip should be subdominant, 
providing a small perturbation over the 1/r potential. This can be checked explicitly by 
perturbatively solving the static equation of motion for ijj, 

V^A + V^/^ (i^i^f) j = j^T . (59) 

For a static point-like source T = — M5'^(x), this reduces to 



4 X 

^' + A^.^^')^.=---. (60) 
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Expanding ip into contributions of different order in powers of the non-linear interaction, 

iP = i>o + i>, + (61) 

witli %jjo = — 4^|^p^ , one finds tliat the first perturbation should satisfy the following 
equatioij^, 

. 4„-4 / M \ ^ 

which is solved by the following expression, 

(2n - 1) . 4n-4 / M \ ^ 2n~3 , , 

V'l = -7 — — r^. 63 

(2n-3) V4vrMpy ^ ^ 

One can estimate the crossover distance r^, for the ip profile as the one for which the 
perturbation theory breaks down, 

M ^ 



This scale is of the same order for any n, up to an weakly n-dependent multiplicative 
factor of order one. The expression for the physical field (f) inside the Vainshtein radius is 



(2n-l)^4n^/ 



M 



{2n - 3) \A7rMp 



(65) 
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